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^ . Abstract 

^ The problem of gravitational fluctuations conflned inside a flnite cutoff at radius r = Tc 

^ ! outside the horizon in a general class of black hole geometries is considered. Consistent 
O . boundary conditions at both the cutoff surface and the horizon are found and the re- 
suiting modes analyzed. For general cutoff the dispersion relation is shown at long 
wavelengths to be that of a linearized Navier-Stokes fluid living on the cutoff surface. A 
cutoff-dependent line-integral formula for the diffusion constant D{rc) is derived. The de- 
pendence on Tc is interpreted as renormalization group (RG) flow in the fluid. Taking the 
^ ■ cutoff to inflnity in an asymptotically AdS context, the formula for D{oo) reproduces as 
I a special case well-known results derived using AdS/CFT. Taking the cutoff to the hori- 
zon, the effective speed of sound goes to inflnity, the fluid becomes incompressible and 
the Navier-Stokes dispersion relation becomes exact. The resulting universal formula for 
the diffusion constant D (horizon) reproduces old results from the membrane paradigm. 
Hence the old membrane paradigm results and new AdS/CFT results are related by RG 
flow. RG flow-invariance of the viscosity to entropy ratio ^ is shown to follow from the 
flrst law of thermodynamics together with isentropy of radial evolution in classical gravity. 
The ratio is expected to run when quantum gravitational corrections are included. 
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1. Introduction 

If you disturb a black hole, the horizon responds like a viscous fluid In par- 

ticular, damped shear waves propagate outward from the disturbance. The fluid viscosity 
was computed to bJ rj = ]^Q^Q in a prescient 1978 thesis by Damour [Q, see also |[5|-[T0[]. 
Dividing this viscosity by the Bekenstein-Hawking horizon entropy density s = yields 
the dimensionless ratio f = The ingredients of the computation are only the basic 
equations of general relativity and (in the computation of s) quantum mechanics. 

A quarter century later, f = ^ reappeared in a beautiful paper by Policastro, 
Starinets and Son (PSS)[0, see also [p!2[ - p8| and the review [|2^. They considered the 



hydrodynamics of supersymmetric A/" = 4 U{N) gauge theory. Stress tensor correlators 
at spatial inflnity were computed in an AdS-black-brane geometry giving 77 = ^N'^T^ and 



Throughout this paper we use G for Newton's constant and set h = ks = 1- 
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s = ^N'^T^. This again yields f = 4^- The ingredients of this computation are much 
more extensive than that of Damour and include string theory, supersymmetric gauge 
theories and AdS/CFT. 

On the face of it, the Damour and PSS computations are very different. However, 
there are strong indications that they are related. Firstly, they both relate a theory of 
gravity to a "dual" fluid theory living in one fewer dimension; without the radial direction. 
Secondly, both approaches lead to the same numerical ratio for ^. An important difference 
is that the Damour calculation is performed at the black hole horizon r = rh, while the PSS 
calculation is performed at spatial infinity r = 00. Both the basic relation between redshift 
and radiu^ and, in the special context of string theory, results from AdS/CFT suggest 
that from the point of view of the fluid theory (which does not have an r coordinate) 
changing r is equivalent to renormalization group (RG) flow. Hence one expects the 
Damour calculation to be related to the PSS calculation by some kind of RG flow into the 

In order to verify this expectation, one 
must flrst deflne what one means by the gravity and fluid theory associated to flnite r. In 
this paper, among other things, we propose a precise deflnition of the flnite r theory and 
show that the expectation is indeed realized. 

The basic idea is to introduce a cutoff surface Ec at some flxed radius r = Tc outside 
the black hole or brane. We then impose ingoing boundary conditions at the horizon and 
flx the induced metric on E^. These boundary conditions do not fully specify a solution. 
The problem is then to identify the remaining internal degrees of freedom and describe 
their dynamics. We solve this to linear order in the internal fluctuations (in appropriate 
expansions) and show they correspond to those of a fluid. A formula is derived for the 
diffusion rate and other hydrodynamic quantities, which generically run as a function of the 
cutoff Tc- We hope that it is possible to extend our approach beyond the next-to-leading 
order considered here, but we defer that problem to future work. 

This reformulation of fluid/gravity duality is the analog - or holographic dual - of 
Wilson's reformulation of quantum fleld theory. Wilson did not insist on an ultraviolet 
completion of quantum fleld theory, and we do not insist on an asymptotically AdS region 

^ Radial transformations are referred to as renormalization already in 

^ These references typically study RG flow by looking at the r-dependence of correlators whose 
boundary conditions are imposed in an asymptotically AdS region. This corresponds to choosing 
a specific UV completion. In this paper we will formulate the problem in a way that does not 
involve such a choice. 



IR. This view is advocated in iT3,B0,P9.31 
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of the geometry. Specifying the couphngs at the Wilsonian cutoff Ky/ is the analog of 
specifying the boundary conditions for the induced metric (and other fields if present) on 
Ec . If we scatter fields at energies below Avk, we needn't know anything about the theory 
at energies above Avk- Similarly if we disturb a black hole by throwing something at it from 
the radius r^, we needn't know anything about the geometry outside Tc-B One advantage 
of the Wilsonian approach is that a much broader class of theories can be discussed. 

In addition to providing new methods of computation and broadening the space of 
applications, several new qualitative insights are gained in this approach. As the cutoff is 
taken to the horizon — > r/^), the relevant geometry is simply Rindler space, and the 
transport coefficients all approach simply computable universal values which are largely 
insensitive to details of the geometry or matter couplings.! In particular — ^ 4^- More- 
over the fluctuations dissipate according to the linearized Navier-Stokes equations, with 
no higher-derivative corrections. Hence our near horizon scaling is the geometric version 
of the low-velocity scaling in which fluids are governed by (as it turns out incompressible) 
Navier-Stokes. At any finite Tc > rh there are infinitely many higher derivative corrections 
to the Navier-Stokes dispersion relation, and one computes only the leading term at long 
wavelengths. The leading dispersion constant in general runs and does not take a uni- 
versal value at radial infinity. If we specialize to an asymptotically AdS black brane and 
take Tc — )■ 00, our computations are all in manifest agreement with the usual AdS/CFT 
definitions of the transport coefficients. All of this supports the picture that Damour was 
computing in the IR of the dual fluid theory while PSS were computing in the UV. The 
extra ingredients required for the PSS calculation are the extra ingredients needed to spec- 
ify a theory all the way up to the UV, while the universality of ^ is a characteristic of the 
IR fixed point and needs only the IR near-horizon Rindler space. 

Interestingly enough, in contrast to the generic transport coefficient, in the classical 
gravity limit the particular ratio ^ typically does not run and equals ^ everywhere. This 
is why the UV and IR PSS and Damour calculations agree for this quantity. Although still 
partially mysterious to us, we show that this RG-flow invariance stems from the first law 
of thermodynamics applied to the radial fiow, together with the fact that, in the classical 

^ In this analogy, large rc corresponds to large Avk, but a precise functional relation will not 
be found herein. 

^ This remains true as long as the linearized gravity fluctuations are governed by the linearized 
Einstein equation, which is not the case with higher derivative corrections. 
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gravity limit, there is no entropy except horizon entropy and the flow is therefore isentropic. 
This will not be the case when quantum corrections are included, as there is then entropy 
in the gas of Hawking radiation as well as entanglement entropy across Ec. Therefore we 
expect ^ to run at the quantum level. 

1.1. Applications 

This work has several distinct applications and motivations. 
• Condensed matter / cold atoms 

One possible application is to the recent exciting work (for reviews see [|36| , |37| , p8| ) on 



trying to solve or understand a variety of condensed matter and cold atomic systems via 
the construction of holographic gravity duals. Most of this has involved asymptotically 
AdS geometries. If this program is to succeed, it cannot depend on the existence of such an 
AdS region (or underlying string theory), as one thing we know for sure about real world 
materials is that they are not conformally invariant in the UV, and certainly do not become 
supersymmetric = 4 Yang-Mills or critical string theory! Rather success of this program 
should depend on the more universal features of the holographic relation between gravity 
and matter systems. The existence of such a universal relation became clear already in the 
70's, but many aspects of its form and origin remain(ed) deeply mysterious. The emergence 
of partially soluble examples within string theory has enabled us to guess much about how 
this more general relation must work. In this paper, we have set up a framework which 
attempts to liberate the essential features of fluid/gravity duality from the confines of the 
string theory/AdS context. We hope that this will enable a broader class of applications, 
and circumvent the often unnecessary discussions of the asymptotic AdS region.! 
• The Navier-Stokes equation 

The Einstein and Navier-Stokes equation are two of the most important differential 
equations in physics and mathematics. Any direct relation between them is extremely inter- 
esting. For example, it was suggested in that the problem of Navier-Stokes turbulence 
might be mapped to a problem in general relativity, with the different scales appearing in 
turbulent phenomena being mapped to different radii in the dual geometry. We hope that 



^ This view has already begun to emerge in various discussions. For example, in a recent talk 
[3S], Polchinski drew a line at the inner (IR) boundary of the AdS region and threw the AdS 
region out because it was irrelevant to the IR physics under consideration. In this paper, we 
attempt to understand how the boundary conditions at Polchinski's line (our Sc) are defined so 
that we need never draw the AdS region in the first place. 
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our approach will facilitate this mapping, as one need not include an AdS region (irrele- 
vant for turbulence), and the forcing term can be supplied by time-dependent boundary- 
conditions on Ec. Of course, for this perspective to be useful, we must understand the 
duality beyond the linearized fluctuations considered in this paper. 
• The gravity /thermo /quantum triangle 

The previous two applications aim to employ gravity to solve problems in other areas 
of physics or math. However it is possible that the ideas involved have more than technical 
utility. Ever since the seminal work of Bekenstein and Hawking, it has been clear that 
there is a deep and fundamental relation between gravity, thermodynamics and quantum 
mechanics, while its detailed form and origin was and is largely mysterious. Much light has 
been shed on this triangle in the specific context of string theory following the microscopic 



description of black holes as a finite temperature two dimensional CFT and its higher- 
dimensional generalization to AdS/CFT It seems likely that the basic triangular 
relationships transcend string theory and AdS/CFT, although lessons from string theory 
are likely useful guides for unraveling the more general picture. 

It is our hope that the attempt here to generalize fluid/gravity duality away from 
the stringy context to its most essential ingredients may be useful in understanding this 
triangle. Much of the current work on fluid/gravity duality attempts to learn about matter 
systems from gravity: we would like to reverse the arrow towards an understanding of 
quantum gravity. In this regard, there may be interesting connections with previous work 
in this direction including in particular ||4^ , p| , |45[ . 



1.2. Outline 

This paper is organized as follows. In section 2 we write down the most general p + 2- 
dimensional general geometry which, on symmetry grounds, could serve as a holographic 
dual to a fluid in p + 1 dimensions. Explicit expressions are given for the asymptotically 
flat and asymptotically AdS black branes with and without charges to serve as illustrative 
examples. 

Section 3 treats the case of an electromagnetic fleld propagating on these geometries 
as a simple warm up. Dirichlet (ingoing) boundary conditions are imposed at the cutoff 
surface r = (horizon r = r^). It is then shown in a long- wavelength expansion that 
the remaining dynamical modes are described by a charge density which evolves according 
to Fick's law in p + 1 flat dimensions. The diffusion constant is given by a line integral 
of certain metric coefficients from rh to Tc, and "runs" as Tc is varied. If the cutoff Tc 
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is taken to r/j, no long wavelength expansion is needed and the Fick law becomes exact. 
Moreover it is shown in this limit that, after carefully normalizing by the divergent local 
Unruh temperature, the diflFusion constant approaches a universal constant determined by 
the properties of Rindler space. 

In section 4 the analysis is adapted to linearized gravity fluctuations. After fixing the 
cutoff and horizon boundary conditions, the vector or shear modes are shown to obey the 
linearized Navier-Stokes equation in a long wavelength expansion, and a running formula 
for the diffusion constant is derived. As the cutoff is taken to the horizon, the linearized 
Navier-Stokes equation becomes exact and the constant is shown to approach the universal 
Damour value. Some special features of the RG flow for the gravitational case are also 
discussed. Tensor modes are shown to have no dynamics in the appropriate limit, while 
there is a dynamical "sound" mode in the scalar sector. Very interestingly, the effective 
speed of sound goes to infinity and hence the sound mode decouples as the cutoff' is taken 
to the horizon. This means the fiuid is becoming incompressible. Specific examples of the 
charged and neutral AdS black branes and asymptotically fiat 5''^-reduced NS5 branes are 
worked out in detail. 

In section 5 we introduce the Brown- York stress tensor tab on the cutoff surface. Prior 
to this point only equations of motion have been used so entropy, viscosity, energy and 
pressure (which depend on the normalization of the action) could not be discussed. We 
show that tab not only is conserved with our boundary conditions but takes the form of a 
fluid stress tensor (to linear order). We compute the thermodynamic quantities in terms 
of the spacetime geometry. 

In section 6 we compute the viscosity to entropy ratio ^ and show that, under rather 
general assumptions, the radial RG evolution equations imply it is cutoff independent and 
equal to ^ for Einstein gravity. It is shown that these radial equations are nothing but - 
in the fluid picture - the flrst law of thermodynamics for isentropic variations in disguise. 
The radial flow is isentropic because in classical gravity there is no entropy outside the 
black hole. It is accordingly suggested that the RG-invariance of ^ will be violated by 
quantum gravity corrections. 

2. Background geometry 

2.1. The general case 

In this paper we are interested in studying the dynamics of fluids in p flat space 
dimensions and one flat time dimension. The holographic dual of such a fluid in its 
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ground state should be ap+2-dimensional spacetime geometry, with isometries generating 
the Euchdean group of p-dimensional rotations/translations plus time translations. The 
corresponding line element can be written in the form 

= -h{r)dT^ + 2dTdr + e^^^'^^dx'dxi, (2.1) 

where the index i = l,...p here and hereafter is raised and lowered with Sij. Lines of 
constant r and varying r are null. We consider the case where there is a horizon rh at 
which 

Hrh) = 0, (2.2) 

and h{r) is positive for r > Vh- Lines of constant r = Vh and varying r are the null 
generators of the future horizon, while those of constant r > th and varying r are timelike 
and accelerated. For convenience we choose the scaling of the spatial a;* coordinate so that 

t{Th) = 0. (2.3) 

A special role will be played by the "cutoff" surfaces Ec of constant r = Vc > r^. The 
induced metric on such a surface is flat p + 1-dimensional Minkowski space 

dspj^-^ = -h (^dr - + e^^dx'dxi, (2.4) 

with T the time coordinate. We will sometimes collectively denote the Minkowskian coor- 
dinates by 

a;"~(a;\r), a = 0,---p. (2.5) 
It is convenient to introduce proper intrinsic coordinates on Sc 



xl = ^/h{rc)T, = e*^^=)a;\ (2.6) 
The advantage of these coordinates is that the induced metric is simply 

dsl+i = riabdx^dx^^, (2.7) 
so that they directly measure proper distances on Eg. Full bulk coordinates will be denoted 

x^" {x\T,r), // = 0, + (2.8) 
We denote by l^^ the normal satisfying 

= [dr + hdr), f = h. (2.9) 
At r = r^,, ^ is null, normal and tangent to the future horizon. 
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2.2. Some special cases 

Here we collect some specific examples which will be used as illustrations in the text, 
of course reduces to flat space for h = t = 0. Another useful way to write flat space 
is in "ingoing Rindler" formS 



h{r) = r, t = 0, 

ds^ = —rdr'^ + 2dTdr + dx^dxi. 

Observers at fixed r > and x^ are then Rindler observers. 
The asymptotically AdSp-(-2 black p-brane solutions are 

2 P+1 



^2.10) 



2 P+ J- 2 

ds'if^ = -T^l 1 ^ Idr^ + 2drdr + -^dxidx'' 



(2.11) 



„2 , ^2 

::2' 

h 

Rindler space (|2.10|) is a limit of the black brane geometry (|2.11J ). To see this define 

, ^ R^{r-rh) 

(P + IK ' (2.12) 

, {p + l)rh 

T = T. 

The horizon is then at r' — i- near which 

ds\B = {-r'dr'^ + 2dr'dr' + dx'dxi) + • (2-13) 

Hence Rindler space is both the near-horizon and — t- oo limits of the black brane. 

If we add a t/(l) gauge field and charge density Q to the black brane, the metric and 
gauge field A are 



e 



^2 I V ^ / ^p+l ^ ^2p ) ' ' 

/X (2-14) 
p - 1 I rP-iy 

Here a = ^^^r^, and we have set the electromagnetic coupling constant to one. 



''' Writing r = 21nt^, r = —t^t , the 2D part of the metric becomes —idt^dt . 
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A specific example with no AdS region is provided by the weU-studied asymptoticaUy 
fiat NS5 brane (we omit the case of general p for brevity). This is a solution of ten- 
dimensional supergravity with a three-form field strength threading an which surrounds 
the brane. This ten-dimensional geometry of course is not of the form (p.l|) . However if 
we Kaluza-Klein reduce to 7 dimensions on the S^, then it does take this form. In the 
7-dimensional Einstein frame and coordinates ( p.l| ), the metric is 



7-2\V5/ 2\ ,,6/5 / r2\l/5 / r2\-V5 



Here y{r) is the solution of 



ry / r2\Vio 

r = y dy' y''/' + • (2-16) 

The right hand side is a hypergeometric function. 

We wish to stress that our approach applies to geometries of the general form ( |2.1| ) and 
is not tied to the above specific examples. Other interesting examples include the proposed 
holographic duals to superconductors , for which the metric cannot in general be found 
analytically, or other cases ( |l47| , ^ to mention a few) which are not asymptotically AdS 
and correspond to systems which are not conformally invariant in the UV. 



3. Electromagnetic warmup 

In this section we warm up to the gravity problem by considering the conceptually 
similar, but mathematically simpler, problem of an electromagnetic field F propagating in 
the geometry (|2.1| ). 

3.1. The setup 

Our first step is to introduce a cutoff surface outside the horizon in the general 
geometry ( |2.1| ) 

Ec : r = Tc > rh, (3.1) 

with coordinates 

x"~(x\r). (3.2) 
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The induced metric on is flat and given by ( |2.4| ) with r = r^. We wish to study the 
dynamics of an electromagnetic field F within the region 

rh <r < Tc. (3.3) 

This requires boundary conditions at both Vc and rh- Since rh is a black hole horizon, we 
impose ingoing boundary conditions there. At Tc, a natural Dirichlet-like choice is to fix 
the components of the field strength tangent to 

Fab{x^,rc) = fabix^) a,b,e = 0,...p. (3.4) 

We view the /s as the parameters defining the cutoff theory, and Ec as the place where 
experiments are set up and measurements made which probe the entire region (|3.3|) below 
the cutoff. Fixing a radius where experiments are performed is dual, in the fluid picture, 
to flxing the scale at which experiments are performed. 

The boundary conditions at Tc and do not uniquely specify a solution of the Maxwell 
equations for F. The problem is to describe the remaining dynamical degrees of freedom. 
We will see that, in the limits we consider, they are described by a single function ^(x") 
which obeys a simple diffusion equation on p+l-dimensional Minkowski space, q can be 
thought of as a charge density on the horizon or, by a simple rescaling, as a charge density 
at the cutoff Qc- We will compute the diffusivity and also see that the data fab specifying 
the boundary conditions at Ec function as a source for the charge density. This enables 
cutoff observers to probe the dynamics of the charge density. 

Let us now turn to the details of this description. 

3.2. Equations and boundary conditions 

We have a bulk gauge fleld with components Fm Fir, Fir, Fij. In terms of these the 
bulk Maxwell's equations may be written 

r : e^^drFrr + d'F^^ = hd'F.r, (3.5) 

i : drF,r + {{2-p)t' -(p')Fri-drFri + h'Fir = -hdrF,r + h{{2-p)t' -(p')Fir + e-'^^d^Fj,, 

(3.6) 

r: drFrr + ipt' + (l}')Frr + e-'''d'F,r = 0, (3.7) 
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where 9* = S'^Wj and for future utility we have aUowed for a position-dependent gauge 
couphng = normahzed so that g{rh) = 1 and (/>(r/i) = 0. In addition we wiU need 
the Bianchi identities 

drFir = drFir — diFm (3-8) 

^rFij d'lFj-j djFj^^. ^3.9) 

We wish to impose ingoing boundary conditions on the gauge field at the future 
horizon. As our coordinates are regular on this horizon, we require F to be regular there:! 

F,,(r,,) = finite. (3.10) 

The other data at the horizon are the horizon current and charge density, defined asi 

= F,,(a;", rh), q{x'') = Frr{x\ r^), (3.11) 

together with Fij{x"',rh). Given the regularity condition, the Maxwell equation ( |3.5| ) at 
the horizon becomes current conservation 



drq + d'j^ = 0. (3.12) 



3.3. Long wavelength expansion 



In this subsection we introduce a non-relativistic long-wavelength expansion which is 
suitable for studying hydrodynamics. 
• Solving the equations 

The general solution of the Maxwell equation cannot be found analytically in a general 
geometry of the form ( |2.1| ). To proceed further we consider a long- wavelength expansion 
parameterized by e — )■ 0. We take temporal and spatial derivatives to have the non- 
relativistic scaling 

~ e^ a, ~ e. (3.13) 



We thank Stephen Green and Robert Wald for pointing out an error in the previous version 
of our paper. 

^ We may also write this in terms of the normal ( |2.9| ) as {q,ji)°' = F°'^lh. 
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The gauge field has the associated expansion 



F,, = e(FO+ei^i + ■■■), 
Frr = e^F^^ + eF^^ + ---), 

(3.14) 

F,,=e2(Fg + e4 + ■■■), 
F,r = e'{F^^ + eFl, + ---). 

We wiU solve for F^^^x"", r) in terms of its value at the horizon r = r/i by integrating the first 
order radial evolution equations (|3.6|) and (|3.7| ) outward to r > r^, and then demanding 



agreement with the boundary conditions ( |3.4| ) when the cutoff r = Tc is reached. At lowest 
order in e these equations are 

zO : h'Fl = -hdrFl + h{{2 - p)t' - (t)')Fl, (3.15) 

rO : drF^, + {pt' + (j)')F^^ + e'^'d'Fl = 0. (3.16) 

There are no non-trivial solutions of ( |3.15D which obey the ingoing boundary condition 
i WM - Therefore 

FO = 0. (3.17) 



The general solution of the second equation (3.16) is then 



F,o^(x^ r) = -e-P*('^)-'^('^)Q°(x"). (3.18) 
The leading order Bianchi identity 

drFl = e-P*(")-'^(")a,Q° (3.19) 
then implies the leading term in Fi^- 

Flix^^r) = - [ dse-P'^'^-^^'^d.q^x") + /° (x"), (3.20) 

J r 

where we have used the boundary conditions at the cutoff to determine the integration 

constants. Evaluated at the horizon r = ( p.20|) gives the Fick-Ohm law in p + 1 
dimensions 
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with diffusivity given by the hne integral 

Df^ir,) = I dse-P'^'^-'^^'\ (3.22) 



rh 



Current conservation (|3.12| ) then imphes 

a^gO = D^^d^q'' - d'fl. (3.23) 

In particular if we choose conducting boundary conditions so that the electric field vanishes 
at the cutoff we find Pick's second law 

a.gO = Df^'d\'. (3.24) 

Taking the Fourier transform of this equation, we see that the charge density propagates 
according to the dispersion relation 

iu = Df^'e. (3.25) 
We still need to solve for F[j{x°', r). The leading term in the Bianchi identity ( |3.9| ) is 

= 0. (3.26) 
The solution of this with the given boundary conditions is simply 

i^S(^",r) = /S(x«). (3.27) 

Hence there are no dynamics associated with FP-. 

To leading order in e, ( p.l7|) , (|3.18|) , ( |3.20|) and ( p.27|) comprise the most general solu- 
tion of the Maxwell equation with the given boundary data. The solution is characterized 
by a single function in p + 1 dimensions, the charge density q^{x°'), which obeys the dis- 
persion relation ( |3.23| ). Away from the scaling limit e = 0, the dispersion relation ( |3.25| ) 
for q has corrections of order k^. 

Other parameterizations of the solution space are possible. One alternative is to 
take the normal component of the electric field at the cutoff, which is not fixed by the 
boundary condition there. This is simply related at leading order to the horizon charge 
density, according to ( p.l8|) , by0 

Q°(x") = e2*('^=)F°^(a;", r,) = e-(p-2W-c)-'^('^=)gO(x"), (3.28) 



We could of course have couched our entire argument in terms of Qc, but the derivation is 
then burdened by extra terms which cancel in the end. 
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and hence obeys the same dispersion relation. This parameterization is more natural from 
our perspective as we want to think of measuring the dynamical modes and dispersion 
constants with experimental devices positioned on the cutoff surface Ec- The charge density 
can be sourced, according to ( p.23|) , by turning on an electric field fir- Hence the dispersion 
constant can be measured by turning on a tangential electric field at the cutoff and watching 
how the normal component decays. We shall henceforth adopt Qc as our basic variable. 
• Normalizing the diffusivity 

As it stands, the value of Df^^ ( |3.25| ) is not very meaningful because it is dimensionful 
and can be set to any value by a change of coordinates. We therefore introduce the proper 
frequency uic and proper momentum kc conjugate to proper time and distance in the cutoff 
hypersurface r = Tc- These are 

u:, = —^^-i^, kl = e-'^''^h' ^i^. (3.29) 
v/%J ' dxl 

Observers at the cutoff equipped with a thermometer measure a non-zero temperature 
Tc.0 To determine Tc, note that very near the horizon r — )■ r/i, observers with worldlines of 
fixed r and are highly accelerated. They therefore detect a Rindler temperature which 
diverges as^l 

drh \/h'{rh) , . X , . 

^« = = ^ ^^^^^ ^ ^ ^ 

for any smooth quantum state. For an equilibrium state, the temperature as a function of 
r is determined by the Tolman relation 

T{r)Vh = Th = constant. (3.31) 

For an asymptotically flat spacetime, /i — 1 at infinity, and Th is the Hawking tempera- 
ture. The relation ( |3.31|) together with the boundary condition (|3.30|) determines the local 
temperature at the cutoff to be 

Tc = T(rc) = -^, Th = ^. (3.32) 



The formulation here naturally leads to a cutoff-dependent temperature. If we want to keep 
the temperature at the cutoff fixed as we change rc, we would have to simultaneously move in the 
space of black-brane solutions. 

We have set h = 1, otherwise it would appear on the right hand side of this equation. 
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The quantity Df^ defined by 



jjEM 

luj, = -^kl (3.33) 



then gives the diffusivity in units of the cutoff temperature. Our final expression for the 
coordinate-invariant, dimensionless diffusivity Df^ is 



For Tc ^ Th, has the universal behavior 

5f ^ = i- + 0(r, - r,) (3.35) 

for any geometry. Unsurprisingly, it follows from ( |2.1U| ) that the leading term is exact 
for all Tc for Rindler space. For an AdS black p-brane with p > 1 and a constant gauge 
coupling {(p = 0) the integral is easily evaluated and yields 

I 1 1 _ ( \P-1 

47r(p-l) l-(^)P+i ^ ^ 

Near the boundary Tc = oo we obtain 

Dg^ = P + \ +o(^-^]. (3.37) 



47r(p — 1) \r 



The leading term agrees with the result of Kovtun and Ritz and Starinets [^. We 
see that the dimensionless diffusivity runs from their result to the universal ^ in ( |3.35| ) 
as the cutoff is taken from infinity to the horizon. 

3.4- Near horizon expansion 

Here we consider the near- horizon dynamics when r ^ r^- There is then no need 
to make the long- wavelength e expansion which required small oj and k. Inspired by the 
diffusion behavior of (|3.25|) , or equivalently (|3.33|) , we let (r, r) — > (t, p) where 



r-p r^rh + Xp. (3.38) 
Taking the limit A — > corresponds to focusing on the near-horizon region. We also have 

h{r) = Xp + 0{X^), e*(^) = 1 + C>(A). (3.39) 
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Since we wish to consider couplings which are regular in r we can write 



g</> = g0o ^ Q^X) (3.40) 
where is a constant. The gauge field has an associated expansion 

F,u = FlS)+F'^^J\ + 0{\'). (3.41) 
At leading order ( |3.6| ) and ( |3.7| ) give 

Flf = 0, d,F,f = 0. (3.42) 

From ( |3.6| ) and (^]^) we find 

9^4'^ = d'F.f, = (p - d,F,f + Ut, (3.43) 



where fa is any function of (t,a;*). The cutoff surface is at Tc = Vh + Xp^ and there we 
impose the boundary c 
(p|) at C(A) leads to 



impose the boundary condition F^j^^ = 0. Identifying F^^^ = e '^°qc together with equation 



dtF.f = p^d^F^^ (3.44) 

which is equivalent to 

dtqc = Pcd'^qc- (3.45) 
Alternatively we can write this in terms of the dispersion relation 

ioj^ = -^kl (3.46) 

with no suhleading corrections for r^ ^ r^. This makes sense since taking the cutoff to 
the horizon should correspond to an infrared limit in which higher derivative corrections 
are scaled away. 

3.5. Suhsummary 

To summarize this section, for large Tc, which should correspond to taking the fiuid 
cutoff to the UV, the charge density obeys a complicated non-linear dispersion relation. 
For small u ~ /c^, Fick's law holds (for any r^) with a dimensionless diffusivity 
expressed as a radial line integral from the horizon to the cutoff. In the infrared limit as 
the cutoff is taken to the horizon, the dispersion relation simplifies to the Fick law and 
takes the universal value 
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4. Gravity 

Now we will adapt this to gravity. The linearized Einstein equations become rather 
complicated in the general background (^3) have been completely analyzed in a series 
of papers whose results we shall use without rederivation. One important feature, demon- 
strated in [|T^ and |T^, is that the linearized gravity fluctuations decouple into so-called 
vector, scalar (or sound) and tensor type, characterized by their transformations under the 
0{p) rotational symmetry. We consider these different types in turn. 

We stress that throughout this section we assume that linearized metric fluctuations 
are governed by equations of motion whose form is given by the linearization of the Einstein 
tensor. corrections or even certain matter couplings could change this. Quite often 
these equations are unchanged, especially for the shear mode. We expect qualitative 
generalizations of our results to hold for any type of couplings. 

4.1. Vector/shear modes 

In this section we consider the vector fluctuations, which turn out to be the most 
interesting for our purposes. These have nonzero hia components which obey 

d'h^r = 0, h,j = F(r)%/i,),. (4.1) 

In terms of these we deflne U(l) gauge flelds 

A' = e-^^hidx"", (4.2) 

and fleld strengths 

F' = dA\ (4.3) 

Here the index i labels the different fleld strengths. Following |T^, by considering dimen- 
sional reduction along the polarization direction of the metric perturbations, the linearized 
Einstein equation for each is precisely that of an abelian gauge fleld with position- 
dependent gauge coupling 

e'^ = e^*. (4.4) 

This reduces the equation for the metric vector fluctuations to p — 1 independent copies of 
the Maxwell equations. We may therefore read the answer off of the solution of the previous 
section. We take ingoing boundary conditions at the horizon and Dirichlet boundary 



17 



conditions on h\{x°'^rc) at the cutoff. This amounts to fixing the induced metric on Ec. 
The analog of the charge density qc is the vector field 

v'ix") = {dr - 2t')hl{x'', r,). (4.5) 

It follows from (|4.1|) that these are divergence free 

d,v' = 0. (4.6) 

The analog of the Fick-Ohm law is the forced linearized Navier-Stokes equation 

drv' = Dcd\' + s' (4.7) 



where the diffusivity is 



rh 



-'''^'^^\djdrh]{x\r,)-d^K{x\r,)). (4.9) 



Here the forcing term Si is 

s' =e 

We wish to define a dimensionless coordinate invariant diffusivityli^ Dc by transform- 
ing to proper coordinates and multiplying by the local temperature T^. This proceeds 
exactly as in the electromagnetic case and the result is 

= ^^'^''[^'^'-'^) r d,^-iP+ms)_ (4.10) 



rh 



When the boundary conditions are chosen so that the forcing term vanishes, the transverse 
velocity fields propagate with the dispersion relationlll 

iuj, = D,kl (4.11) 



The quantity Dc, which is the kinematic viscosity times the temperature, differs from the 
dynamic viscosity 77 by a factor involving the energy density, temperature and pressure. It is a 
more basic quantity from our perspective in that it is directly related to the measured decay rate 
of a perturbation. At this point we cannot define r], because we have not determined the energy 
density or pressure. This requires a bit more machinery - the Brown- York stress tensor - and will 
be worked out in section 5. 

We have set h = c = 1, and the Newton constant G does not enter our calculations. In 
more general conventions iuJc = hc^Dckc- The h appears in this classical calculation because the 
diffusivity is expressed in terms of the local temperature which is itself proportional to h. 
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The dimensionless diffusivity behaves universaUy in the infrared as the cutoff is taken 
to the horizonlll 



^ for Tc rh- 



(4.12) 



This agrees with the result obtained three decades ago in [^. It apphes for any geometry 
of the form ( p.l|) . It is reaUy a property of the hnearized Einstein equation in Rindler 
space, which is the only relevant region for the calculation in the Vc — ?■ rh limit. 

4.2. A D-theorem and other special properties of Dc 

So far the discussion of the electromagnetic and gravitational diffusivity have been 
exactly parallel. However some special features arise in the gravitational case. This first 
is that the integrand of ( [4.8| ) turns out to obey 



he 


-{p+2)t- 




-2t'h_ 




_h' 



{h' -2t'hy 



(4.13) 



Here T^;^ is the bulk matter stress tensor and ( is any null vector tangent to the brane with 
time component Hence if there is no matter or if T„,^C^C^ = 0, the integrand is 



a total derivative. Dc is then given by the simple expression 

4n h'{r,)-2t'{r,)h{r,) 



(4.14) 



The fact that the expression for the diffusivity can be integrated stems from the fact 
that the shear modes are pure gauge at zero momentum (we generalize here [^). hir is 
nonzero for /c* = a; = 0, and we can derive its radial dependence from the zero momentum 
Einstein equation. However we can also solve this equation with a gauge transformation 
of the form 



Sx' - eV + Be' 



-2t 



dr, 



6r = -Bcix'' 



(4.15) 



with S, e* arbitrary constants. The second term in Sx"^ is added to preserve our gauge 
condition hr^ = 0. The nonzero component is 



h^r = {Bh + e^*) Ci 



(4.16) 



The results of [1^,2C,24,2^,28| suggest the universal value will be modified by higher- derivative 
gravity corrections. 
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In order to preserve the boundary condition hir{rc) = 0, we must take 
resulting in 

hUr) = e.ie"'^^^ - ^/'^''^) ^ h{r)e,- (4-18) 

Of course we already found this solution, which is not pure gauge at non-zero-momentum, 
using the Einstein equation. The present derivation has the distinct advantage that it 
follows from gauge invariance and therefore is completely independent of the prescribed 
dynamics. Note that we can also find matter perturbations compatible with this gravity 
solution by the same method. 

Dc can be expressed as a ratio of the coefficients of drhir and drhij. The latter 
vanishes at zero momentum so we need to work to next order in e expansion. We promote 
e* — )■ e*(x°'), with die'^{x°') = 0. By symmetry hij = e^*/2(r)(9iej +djei) for some function 
/2(r), determined from the ij components of the linearized Einstein equation 

SG,j = SnGSTij. (4.19) 

Again by symmetry the variation of Tij with the matter fields is of the form tmidiCj +djei) 
and ( [4.19| ) reduces to a single equation which becomes 

e-(P-2)*a^[-e(P-2)*/i + e^^hdrf2] = -IGnGt"^ (4.20) 

in the long wavelength limit. In many cases (all that we have studied) t^, = and hence 

h{r)drf2ir) = e-2*W/i(r) - h{rH)e-P'^'^\ (4.21) 



This is then enough information to reproduce (|4.14|) . So the universality of ( ^1.14|) is at 



least in part a consequence of zero-momentum gauge invariance. 

Interestingly, if we plug in the neutral AdS black brane metric ( p.ll|) we find the 
surprising result 

D^ = l- (4.22) 

for any value of Tc. Hence in this particular case Dc (unlike Df^ for the same geometry) 
does not run. In general one expects Dc and as well as all transport coefficients to run. 
Indeed for the charged black brane ( |2.14| ) we will exhibit a rather nontrivial RG flow in 
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section 4.5 below. The constant value appearing in (|4.22 ) is directly related to the fact 



that ^ does not run, which will be discussed in section 6 below. However for now we briefly 
note the relation. The viscosity is related to Dc by the relation rj = '^'^(^+^) where £ {V) 
is the energy density (pressure). In the absence of chemical potentials, the entropy density 
is determined from TcS = £ + V, implying Dc = ^. When there are chemical potentials 
such as for the charged black brane this simple formula no longer holds, and Dc runs. 
We now show that, in a wide range of circumstances, Dc decreases with increasing r^. 



i.e. it increases under RG flow into the infrared. We follow in spirit the A-theorem of |^ 
Define the quantity 

H{rc) = e^P-^^'^'-^^h\rc)dr(^^)=e-^'h-{h'-2ht')eP' P dse-^^+^W^). (4.23) 
Since h{rh) = 0, we have H{rh) = 0. The gradient of H then obeys 

drH{r) = -IGnGT^^C^Ce^^^'''^ j rfse-(P+2)t(.)_ ^4 34) 
This relation employs the Einstein equation 

IQ-kGT^X^C = '^G^^CC = h" + {p- 2)h't' - 2t"h - 2pt'^h. (4.25) 

The null energy condition implies T^j^C^C^ ^ O5 ^ind hence that drH < 0. Since H{rh) = 
we conclude that H{r) < for r > r^- It then immediately follows from the definition 
O) that 

dr^Dc < 0. (4.26) 

We expect that the general approach of this proof will be applicable to a wide range of 
situations. However, we wish to note that the precise result (|4.26|) may also be obtained 
from (f4.14 ) simply by differentiating with respect to Vc- 



4.3. Tensor modes 

The tensor modes of the metric have nonzero components hij with /i^* = dih'j = 0. 
The equations governing their fluctuations appear in |]T^, and are equivalent to a (p + 2)- 
dimensional scalar Laplacian for each of the {p — 2){p + l)/2 tensor components. The 
analysis of these modes in our setup is very similar to the one given in section 3 for the 
Fij modes of the electromagnetic fluctuations, and will not be spelled out here. The 
conclusion is that, in either the near- horizon or long- wavelength expansions ( p.l3| ), their 



values are everywhere flxed by the boundary conditions on Ec and the horizon. There are 
no dynamical modes in these expansions. 
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4.4- Scalar/sound modes 

In this subsection we consider the scalar or sound mode. The equations of motion 
are somewhat comphcated and we restrict ourselves here to the case of the AdS black 
threebrane so that p = 3, although we expect the more general case to be similar. For 
p = 3 the equations have been fully analyzed in |]T^. They found that the metric fluc- 
tuations are determined by a certain linear combination of the components, denoted Z2, 
obeying a second order radial differential equation which they express (equation 4.35) in 
Schwarzschild-like coordinates. For our purposes it is more convenient to work in the ad- 
vanced coordinates ( |2.11 ), because these are smooth at the future horizon and the ingoing 
boundary condition is simply regularity at r = r/^. Z2, which is a Fourier transform with 
respect to Schwarzschild time, is traded in these coordinates for the Fourier transform with 
respect to r denoted here by X. X is related to Z2 via a factor of e'^ J ^ : 



X 



(4.27) 



Rewritten in terms of X and the coordinate r the equation for the sound mode, equation 
(4.35) of reference [T7\, is 



(r^ - rt) dlX + (r^ - 



5r 



+ 



) rr^/c^ + 3r^ (— /c^ -(- a;^) 



+ 



9 2 



4(r2 -|- r 



drX 



X. 



(4.28) 

Let us now analyze this equation in the long-wavelength expansion (|3.13| ) with uj ~ 
/c2 ~ e2 — )■ 0. Note that there are no poles appearing in ( 4.28 ) for small uj ~ /c2, so we can 



safely take e — )■ for all r. One flnds the equation for the leading term in the e-expansion 
oiX 



r{r 



4 _ r(3r4 - r^) 



drX^ + 



r2(3r'* — r^)(r'^ — r^) 



(4.29) 



which does not depend on a; or /c. This has a unique solution, up to an overall scale, which 
is non-singular dX r = r^. The scale is then fixed by the boundary condition at the cutoff 
Tc- Unsurprisingly we learn that there are no dynamical degrees of freedom in the sound 
mode in the expansion ( |3.13| ). 
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It is interesting, but outside the scope of this paper, to consider an alternate scahng 
hmit in which a; ~ /c. In this limit nontrivial sound modes may appear with some fixed 



velocity Vg- In the limit (|3.13| ), all fixed velocities including Vg are sent to infinity. 

Now let us now consider the near horizon expansion. In this expansion we take — )■ 
without any preassumed relation between ui and k. In this case it is not so simple to take 
fc ^ fh, as a quick inspection of ( [4.28| ) indicates there may be poles at 

u; = ±^k. (4.30) 



Let us first consider the case where ui does not take the value ( [4.30| ). Then we can safely 



take Tc — )■ Vh, and as in our long- wavelength expansion above there are no degrees of 
freedom. Hence the only possibility for dynamical modes are those that obey ([4.3 



Equation ( |4.30[ ) can be written in terms of the coordinate-invariant proper quantities and 
local temperature as 

This is a dispersion relation for a sound mode with velocity of sound 



3 Th 



(4.32) 



However note that as ^ r^, and hence the velocity of sound goes to infinity. Hence, 
in an expansion in r^ — r^, no sound modes appear and the fluid at the horizon is incom- 
pressible. 

This again is consistent with the expectation that the limit — ?■ r/j is a nonrelativistic, 
low energy limit. It is low energy because of the high redshifts, and non-relativistic because 
the degeneracy of the induced metric on appears only in the temporal and not the spatial 
components. In its most general form, the Navier-Stokes equation for a fluid contains both 
sound and shear modes. However one may take a further limit of these equations of the 
form (|3.13|) in which velocities scale as f * — )■ ef* — )■ (see for a nice discussion). In this 



limit the sound velocity goes to inflnity. The fluid retains its nonlinear interactions and is 
described by the incompressible Navier-Stokes equation. The near- horizon limit resembles 
a bulk version of this limit. 

This leads us to the interesting conclusion that the fluid which lives at the horizon is, 
at linear order, universally given by an incompressible fluid with dimensionless diffusivity 



Dc — 
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4-5. Charged black brane 



In this section we consider the case of the charged black brane geometry ( |2.14| ) as 



a somewhat more non-trivial illustration of our approach. If symmetry allows, a metric 
perturbation can source a matter perturbation already at the linear level. For the charged 
black brane there is a matter perturbation of the general form 

5 A = ajdx^. (4.33) 

We can solve the coupled equations by transforming them back to the previous case. First, 
let us define a shift of the "gravitational" field strength appearing in ( [4.2| ) and ( [4 .31 ) by 

F^^ = F^^ - 16nGa^e-^^A'^{r), (4.34) 



with Ar given in equation (|2.14|) . (Note that, in a notational clash, neither F nor F here 



is the field strength of A\ Instead, F refers to the "gravitational" field strength as defined 
in ( |4.2D and (^73|).) F-i^ then obeys exactly the same equations ( P75| ) through ( |3.7| ) obeyed 
by F-/.^ in the neutral case. 

Consequently, following the logic in Section 4.3, we find at lowest order, similarly to 

(EH), 

F^^ = -e-(P+2)*i)J(x") (4.35) 

where 

{)^(x") = F/,(a:^r^). (4.36) 
Note that the effective "charge" in ( |4.36| ) is modified from the corresponding equation 

Now we must examine the two remaining equations: the lowest order Bianchi identity 

= (4.37) 
and the Maxwell equation for the gauge perturbation 

Fi^ = ^dr (^e^P-^'^'hdra'^ . (4.38) 

Note that we have written these two equations in terms of F^^; by doing so we can now 
quickly see that the right hand side of (|4.37|) is just a total derivative. Thus we find, 
similarly to (|3.43| ) for r = r^, 



i r dr dr U^-'^^'hd^dra^) + /,,(x«). (4.39) 

Q Jrh ^ ' 



24 



Before we can use ( [4 .391) to find the diffusivity we must first find . To do so, we 
simply plug ( ^4.35| ) and ( J4.38[ ) into ( [4 .341 ), additionally imposing the boundary condition 
a^ifc) = 0. We thus find 



(r) 



IGttGQ 



.pt(r)-pt(r.) h'{r)-2t'{r)h{r) 



h'irr 



2t'{r,)h{r, 



(4.40) 



Again we can find a dimensionless coordinate invariant diffusivity which is 



d9^ 



-pt{r 



An h'{r,)-2t'{r,)h{r,) 



(4.41) 



p-1 



Note that the first line of (|4.41|) is identical to ( f4.14|) , corroborating the universality of 
( [4.14|) . However when we plug in the metric coefficients for the charged black brane to 
get the second line of ( |4.41| ), we see that, unlike the case of the neutral black brane, the 
diffusivity is no longer constant. 



4-6. Asymptotically flat black p-brane 

In our next example we consider the asymptotically fiat S'"^-compactified black NS5 
solution (p.l5| ). The 7-dimensional effective action for the rotationally invariant modes is 
Einstein gravity plus scalars, hence according to section 4.2, dr^Dc = 0, and Dc = ^ = 
at any scale. Of course this can be reproduced by direct calculation. The local temperature 
is a non-trivial function of radial position: 

-1/2 ^ -1/10 / „,2\ -1/2 

(4.42) 



1 



3/5 



1 + - 



L2 

y'h 



L2 

1 + — 

Vc 



Vc 



General formulae from the following section for the energy density plus pressure give 



£ + V 



18/5 



1 + 



vl 



-3/5 



1 - 



vl 

Vc 



-1/2 



(4.43) 



and entropy density 



-5t(r,) 

4G 



vl 

AGyl 



yi 



-1/2 



1 + - 



L2 



1/2 



(4.44) 



We note that £ + V = TcS. The L — t- oo throat limit of these expressions should describe 
the thermodynamics of the quantum theory on the NS5-brane, but we will not further 
pursue this here. 
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5. Brown- York stress tensor 



Up until this point, all of our results have followed from equations of motion and 
boundary conditions. These are sufficient to determine diffusion rates. However they 
are not sufficient to determine quantities such as the energy, entropy or viscosity which 
depend on the normalization of the action. In this section we introduce the Brown- York 
stress tensor tab whose overall coefficient depends on the action. We propose that, with 
our boundary condition, tab evaluated at the cutoff is the stress tensor of the fluid. As 
we will see, the introduction of tab makes the calculations and analysis considerably more 
concise, as well as enabling us to compute more quantities. 

The Brown- York stress tensor associated to a hypersurface E with unit normal n is 
defined by 

tab = {jabK - Kab + Cjab) , (5.1) 

where 

7ab = 9ab - naUb (5.2) 

is the induced metric on E, the extrinsic curvature is 

Kab = 2^n'~fab 

= I {n^dc-iab + daU'^'Jcb + dbTl'^'^ac) (5.3) 
= \ {^n9ab - riarfCnQcb " nbU^'Cngac) , 

and Cri is the Lie derivative along n. tab is ambiguous up to a constant multiple of the 
induced metric which we denote here C. This is equivalent to the ambiguity in shifting 
the pressure of a fluid by a constant. 

The Einstein equation G^^, = StiGT^^, implies that tab obeys the conservation law 

= n^T^b, (5.4) 

where a, b are projected to E, V is the 7-covariant derivative and T^i, is bulk stress tensor. 
We are interested in the case for which E = Ec is the cutoff hypersurface and will al- 
ways impose "energy-conserving" boundary conditions so that the right hand side of ( |5.4| ) 
vanishes. We then have 

VXfcls. = 0. (5.5) 
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Note that if we were to allow fluctuations in both the induced metric and extrinsic 
curvature, the linearization of ( |5.5| ) about the background would be 



[{6V'')C + V"'{StaB)]^ =0. 



(5.6) 



This equation implies that the linearized Brown- York stress-tensor is not in general con- 
served in the background geometry, precluding a simple interpretation as the linearized 
fluid stress tensor. However with our boundary conditions the induced metric on Ec is 
held flxed and all is well. 

5.1. Background stress-energy 

Now we evaluate tab for our general metric ( |2.1| ) at the cutoff hypersurface Ec. The 
unit normal is 

n = \fhdr + -^dr (5.7) 

and 



K^^dxf'dx'' = Vh 



dr\ 



h { dr — — j + e^^dxidx 



(5.8) 
(5.9) 



In intrinsic coordinates to as deflned in (|2.6|) 



Kabdx^dx 



h' 



{dx'^^f + Vht'dxcidxi, 



2Vh 

'fabdx^'dx^ = rjabdx'^dx\ 
The leading order Brown- York stress tensor is 



(5.10) 
(5.11) 



tl^^dx'^dx'' = ^[-pt'idxlf + {{p-l)t' + ^)dx,,dx'^) +C'7jabdx''Jxl, (5.12) 
ottCt Zri 



where all r-dependent quantities are evaluated at r = r^. This is the stress tensor of a fluid 
at rest with constant pressure V and energy density S. The constant part of the difference 
S — V depends on the choice of constant C: the behavior of the linearized fluid depends 
only on the sum 

^/h / h' \ 

(5.13) 



£ + V= \ —-t' 



SttG V 2h 

Note that for rc ^ £ + V ^ TcS, where s is the entropy density of the horizon. 
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5.2. Perturbations 

Now we consider a small perturbation h. It is convenient to choose the gauge 



hrr — hrr — hri = 0, (5-14) 

SO that the metric retains the general form ( |2.1| ). Moreover, our boundary condition implies 

hatirc) = 0. (5.15) 

It then follows that the leading correction to the extrinsic curvature in the long wavelength 
limit simplifies to 2i^r^^ = \/lidrhab, and 

\/Ti 

tlbi^",r,) = -^{-drhab+-fabl'"'drKd). (5.16) 



For the shear mode, only hij- and hij are nonzero. Moreover, as shown in |T^ the Einstein 
equation implies they are related by Fd(^jhi^j. = hij for some function F(r), as in equation 
( [4.1| ) previously. Conservation of the Brown- York stress tensor implies that F{rc) = ^ 
on Tic and 

drhij = ^drd(^hj)^. (5.17) 



Defining the velocity field 



Vi = TZT^ — :^drhir, (5.18) 



The linearized stress tensor can be written in intrinsic coordinates 

tl^dx^cdxl = 2{S + V)vidxldTc + r]d^Vjdxldxi (5.19) 
with viscosity given by 

(5.20) 



Using (|3^ , (plD and ( |5J3| ), this becomes 



g-pt(rc) 

¥r.) = -j^, (5.21) 

in agreement with Q for Tc — ?■ rh where t — > 0. 

We see from the above that the stress tensor for the linearized shear mode indeed 
takes the form of an incompressible fluid. 
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6. RG flow, the first law of thermodynamics and ^ universahty 

In this section we want to show that the first law of thermodynamics together with 
isentropy of the radial flows is equivalent to a radial component of the Einstein equation, 
and moreover imply that ^ does not run. A general type of equivalence between the 
Einstein equation and the first law has been demonstrated by Jacobson [Q, see also [|5C 



We suspect our equivalence is related to this - as well as the recent work - but we defer 
this issue to future consideration. 

For the present purposes, it is convenient to consider a quotient of the general geometry 
T7\\) under shifts of 



x'r^x'+n\ n'eZ. (6.1) 
This turns the spatial Rp in the p-torus with r-dependent volume 

Vp = eP\ (6.2) 

The total entropy S = sVp as a function of the total energy E = £Vp, pressure V, charge 
Q, chemical potential /i, inverse temperature (3 = and volume Vp are related by 

S = PE + pVVp - P^Q. (6.3) 

In general there will be one nQ term for each chemical potential present; in our case we 
consider only the charged brane given in (|T|), for which n = A^/Vh. 

Now let us consider the first law of thermodynamics for radial variations under the 
assumption that the variation is isentropic. This reads 

= drS = drP {E + VVp -^Q)+P {drE + drVVp + VdrVp - dr^iQ) . (6.4) 

Using expressions ( p.32|) , ( ^.131) and ( |6.2|) for the thermodynamic quantitie^l one finds 
( |6^) becomes 



= J^TTn^^^" + ^P- 2)^'^' - 2^"/^ - ^Pt"h) - ^T.^CC- (6.5) 

ibTVLrlH J- H 

Here we have also used VpT^j^C^^C,^ = Q^e"^* for the charged brane. Comparing with 
( [4.25|) , we see that the right hand side is exactly a component of the Einstein equations. 



We are considering here only the thermodynamic quantities associated to the background. It 
would be interesting, but beyond the scope of this paper, to understand if this reasoning applies 
at quadratic order in the fluctuations, where shear dissipation produces entropy. 
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Therefore, isentropy of the RG flow imphes a radial Einstein equation. Of course this can 
be turned around to state that the radial Einstein equation implies the isentropy of RG 
flow. 

Now if S = constant, then 

^ = IG • 

where the overall multiplicative factor is fixed by demanding the Bekenstein-Hawking law 
s = ^ at the horizon. Combining with ( ^.21| ), we deduce that on Ec 

7^ = r («-^' 

s{rc) 4:71 

for any value of rc.0 We have already seen in our formalism that (6.7) for the special 



value Tc = r/i is a universal feature of Rindler space. Now we see that under quite general 
assumptions, it does not change under RG flow, and so the value will also apply to 
Tc = oo, in agreement with findings in [ pT| , p^ . So far we have not asked the question: why 
should radial evolution be isentropic? A physical answer from the gravity side is that the 
only entropy associated with a classical solution is the horizon of the black hole. Therefore 
the total entropy inside any hypersurface outside the horizon should be independent of 
radius. 

This also gives a clue as to when the relation (|6.7|) might be violated fS^]. When 
quantum corrections are included on the gravity side, the entropy will generically depend 
on radius. For example there might be a thermal gas of Hawking radiation surrounding 
the black hole or entanglement entropy across Ec. These are suppressed by a factor of h 
relative to the horizon entropy. We then see no reason to expect that ( |6.7| ) should survive 
such quantum corrections. The universality of (^]^) is presumably a statement about the 
classical gravity limit. We note of course that the classical gravity limit is in general not 
the same as the classical limit of the quantum theory underlying the fluid. 
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Closely related observations were made in a different formalism by Iqbal and Liu [pi|. 
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